Interferometric description of optical nanomaterials 
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We introduce a simple analytical model that describes the interaction of light with optical nano- 
materials in terms of interfering optical plane waves. In this model, a nanomaterial is considered to 
consist of planar arrays of nanoparticles. In the analysis, each such array reduces to a thin homoge- 
neous sheet that has one transmission and in general two reflection coefficients. These coefficients 
are found to be equal to those of an isolated planar array of the nanoparticles. The presented theory 
enables fast analytical calculation of the light transmission and reflection by nanomaterial slabs of 
arbitrary sizes. It is also shown to accurately describe optically anisotropic nanomaterials that in 
addition exhibit strong spatial dispersion, such as bifacial nanomaterials. 



Optical nanomaterials are man-made materials com- 
posed of densely packed subwavelength-size nanoparticles 
appearing like artificial atoms to light. While the optical 
response of each individual nanoparticle can be revealed 
using, e.g., the electromagnetic multipole expansion^, the 
description of the macroscopic optical response of a real 
three-dimensional nanomaterial still remains a challenge. 
This description is complicated by non-trivial interac- 
tions between the nanoparticles, including evanescent- 
wave coupling between them. In this work, we however 
propose a way to calculate the transmission and reflec- 
tion characteristics of essentially arbitrary nanomaterial 
slabs, including anisotropic and spatially dispersive ones, 
without resorting to evanescent waves. 

So far, both numerical 2 5 and experimental 6 8 tech- 
niques have been used to obtain the transmission and re- 
flection coefficients for nanomaterial slabs with a rather 
limited number of nanostructured layers. The coeffi- 
cients obtained for such slabs do not accurately describe 
the properties of a bulk nanomaterial. One approach 
to extract the bulk properties would be to successively 
increase the number of layers and see if the optical char- 
acteristics converge 2 -"—. However, the understanding of 
the real physics that determines the final transmission 
properties is lost when using this procedure. In the de- 
scription proposed in this work, the properties of a bulk 
nanomaterial slab are directly linked to the properties of 
a single layer of the material. 

Recently, in connection with metamaterials, several re- 
trieval procedures have been introduced to obtain effec- 
tive wave parameters, such as electric permittivity and 
magnetic permeability, from the reflection and transmis- 
sion coefficients of a nanomaterial slab^iS. These re- 
trieval procedures rely upon the Fresnel coefficients which 
are derived for dipolar media. However, for the class of 
bifacial nanomaterials in which the electric quadrupole 
excitations are present 11 , the classical electromagnetic 
boundary conditions do not hol d 12 i 13 . Consequently, nei- 
ther the Fresnel coefficients nor the retrieval procedures 
can be applied to these nanomaterials. The development 
of an adequate theory for the description of highly spa- 
tially dispersive nanomaterials, such as bifacial nanoma- 
terials, would be of great practical importance, e.g., for 
solar cell applications. 



In this work, we re-examine the propagation of light 
through an arbitrary slab of bifacial optical nanomate- 
rial. We find that both the transmission and reflection by 
the slab can be surprisingly accurately described in terms 
of propagating optical plane waves only, which dramati- 
cally simplifies the description. Using rigorous numerical 
calculations we verify the accuracy of our approach and 
find that in many cases, the role of the evanescent waves 
for the energy transfer between adjacent nanomaterial 
layers can indeed be neglected. 

Consider a nanomaterial slab that is created by stack- 
ing two-dimensional arrays of nanoscatterers in a trans- 
parent dielectric medium (see Fig.[]Ji). When illuminated 
by an optical plane wave, the first layer of nanoscatterers 
transmits a certain portion of the incident field. A part 
of this transmitted field is then reflected back by the sec- 
ond layer and the rest is partially absorbed and partially 
transmitted further to the third layer. Provided that 
the nanomaterial periodicity is much smaller than the 
illumination wavelength, no diffracted orders will occur 
and the propagating field between the neighboring layers 
will consist of two optical plane waves. Neglecting pos- 
sible evanescent- wave couplings between the layers, we 
describe each layer as an infinitesimally thin sheet sur- 
rounded by the host dielectric (see Fig. [TJd). The total 
transmission and reflection of the slab is then described 
in terms of the transmission and reflection coefficients of 
the individual sheets in a way resembling the description 
of Fabry-Perot interferometers. 

Our approach is as follows. For a single isolated layer of 
a nanomaterial, we first numerically calculate the trans- 
mission and reflection coefficients and assign them to an 
equivalent infinitesimally thin sheet in the middle of the 
unit cells. These coefficients depend on both the angle 
of incidence and the polarization of the incident field. 
In addition, for bifacial nanoscatterer arrays, the reflec- 
tion coefficient changes if the illumination direction is 
reverse d 11 ! 14 . Therefore, for the wave propagating to the 
right (left) within the nanomaterial, we use the reflection 
coefficient pL (pr) to describe the reflection from the left 
(right) side of each layer. Optical reciprocity ensures that 
the transmission coefficient r must be the same for both 
waves. When calculated in this way, the parameters p^, 
Pr and r automatically include the near- and far-field 
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FIG. 1. A nanomaterial slab (a) that is composed by stacking 
layers of arbitrary scatterers, S, is described as (b) an array 
of infinitesimally thin sheets. The transmission and reflection 
coefficients (r, pL and pr) for each of these sheets are taken to 
be equal to those of a single isolated layer of the nanomaterial. 
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FIG. 2. Light propagation through a nanomaterial described 
by an array of infinitesimally thin sheets. Between each pair 
of such sheets there are two counter-propagating plane waves 
with transverse field components Uj and Uj. The unit cells 
of the nanomaterial are shown with dashed lines. 



coupling between the scatterers within the layer in ques- 
tion. This coupling makes the extent of the evanescent 
wave along the layer's normal shorter for denser packing 
of the scatterers. 

Let a plane wave with a transverse field amplitude Uo 
and a wave vector k[ n = xk x + zk z be incident on a 
nanomaterial slab that consists of N layers of thickness 
A and has its surface normal along the z axis (see Fig. [2]). 
The wave can be either TE- or TM-polarized. Treating 
each nanoscatterer layer as an infinitesimally thin sheet, 
we consider the counter-propagating waves between the 
sheets to have wave vectors k± = xk x ± zk Z: because 
the material between the sheets is considered to be the 
same as outside the slab. The transverse fields Uj and 
Uj after each sheet j in Fig. [2] must satisfy the following 
equations 

U j = fU j . 1 +g R U' j , (1) 



U'j = g L Uj + fU' j+1 , (2) 

where / = rexp(^A), g L = p L exp(ik z A) and qr = 
PRex.p(ik z A). Using Eqs. (pQ) and (j2j), we derive separate 
relations for the forward and backward propagating fields 

Uj+i + Uj-! - aUj = 0, (3) 
U' j+1 + U 3 _ x - aUj = 0, (4) 

where a = f + - gLQR). With the help of Eqs. (pQ) 

and (|3j) and the fact that U N = 0, we obtain the trans- 
mission coefficient of the slab to be 

Uo G N (a) - /Gn-^oY {> 

Here we have introduced the G-polynomial that is calcu- 
lated as 

G o (a)=0, (6) 
G 1 (a) = l, (7) 
Gj(a) = aGj-i(a) — Gj-2(®). (8) 

Similarly, using Eqs. (pQ), (|2]) and (UJ) we derive the reflec- 
tion coefficient 

r = ^ = gLG N (a)f-H. (9) 

Equations ([5]) and (|9j) enable direct calculation of the 
transmission and reflection coefficients of an arbitrarily 
thick nanomaterial in terms of the transmission and re- 
flection coefficients of an isolated monolayer of the nano- 
material. For N = 1, Eqs. ((5]) and ((9]) correctly yield 
t = r ex.p(ik z A) and r = pL ex.p(ik z A) and for N = 2 the 
well-known results for a Fabry-Perot etalon are obtained. 

In order to demonstrate the applicability of our the- 
ory, we compare it with rigorous numerical calculations. 
This is done by selecting some non-trivial nanoscatterers, 
such as nanoshells (Fig. [3^l), nanorings (Fig.H^i) and nan- 
odimers (Fig.|5k). These nanoscatterers are considered to 
compose stacks of two-dimensional periodic arrays that 
are embedded in a dielectric host medium of refractive in- 
dex 1.5. The transmission and reflection coefficients for 
a single array are computed numerically using the com- 
puter software COMSOL Multiphysics. The obtained co- 
efficients are then used in Eqs. (j5]) and (j9j) to acquire the 
transmission and reflection coefficients for several layers. 
These coefficients are compared with the results of direct 
numerical calculations of the whole stack with COMSOL. 
For these calculations we choose TM-polarized illumina- 
tion with = 45° and a slab consisting of N = 5 layers. 
This choice is general enough for demonstrating the ap- 
plicability of the model. 

We first consider a nanomaterial with an isotropic unit 
cell containing a silver nanoshell as depicted in Fig. [3£i. 
The nanoshells have an outer radius R = 35 nm and 
a thickness h = 7 nm. They form a cubic lattice with 
period A = 130 nm. The optical properties of such 
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FIG. 3. (a) Geometry of the nanoshell [R — 35 nm, h — 7 nm, 
A = 130 nm]. (b) Transmission and reflection spectra of a 
five-layer thick slab for TM-polarized incident light with the 
angle of incidence of 45°. The numerically calculated spectra 
(solid lines) are shown along with the the analytical results 
(stars) obtained from Eqs. J5} and ©. (c) Reflection spectra 
for increasing number of layers as obtained from Eqs. (|5]) and 



nanoshells are well investigate d 15 ^ 6 and similar struc- 
tures can be relatively readily fabricate d 17 ^ 8 . For the cal- 
culations, the optical characteristics of silver were taken 
from Ref. [HI. The calculated transmission and reflection 
spectra in the wavelength range from 500 to 1000 nm 
for the nanoshell slab are shown in Fig. [3)3. The spectra 
obtained by using Eqs. ((5]) and (j9j) are in a good agree- 
ment with the direct numerical calculations, indicating 
that the plane-wave description of the light propagation 
is appropriate and that the evanescent- wave coupling be- 
tween the layers is weak. While direct numerical calcula- 
tions are always limited by the present computational re- 
sources, the introduced theory allows us to calculate the 
response of an arbitrarily thick slab. In Fig. [3b the reflec- 
tion spectrum for a slab of 1000 nanoshell layers is shown 
by the red curve. For this particular case of nanoshells, 
the spectrum of 1000 layers is close to the spectrum of 
the 5-layer slab and it is already indistinguishable from 
the spectrum of an infinitely thick nanomaterial. 

Next, we introduce an anisotropic (uniaxial) unit cell 
containing a silver nanoring as depicted in Fig. [4^. The 
ring has an outer radius R\ — 20 nm, inner radius 



FIG. 4. (a) Geometry of the nanoring [R± = 20 nm, R2 = 
10 nm, h = 10 nm, A = 50 nm]. (b) and (c) are as in Fig. [3] 



R2 = 10 nm and thickness h = 10 nm. The rings form a 
cubic lattice with period A = 50 nm, such that each layer 
is aligned with the xy-plane. The calculated transmission 
and reflection spectra for the nanoring slab are shown in 
Fig. [4)3. The theory yields excellent agreement with di- 
rect numerical calculations also for these nanoscatterers. 
The reflection for 1000 layers, depicted by the red curve 
in Fig. [3fc, shows that the bulk nanomaterial behaves 
quite differently from a single layer due to the interlayer 
interaction. One can notice that if the number of layers 
is large, the nanomaterial acts as a spectrally selective 
broad-band reflector with a nearly flat-top spectral pro- 
file. 

As a final example, we consider a bifacial nanomate- 
rial slab that exhibits strong spatial dispersion. To our 
knowledge such nanomaterials have not been previously 
described in terms of their reflection and transmission 
characteristics. The unit cells of the material contain 
asymmetric silver nanodimers (see Fig. [5ji). These nan- 
odimers have been shown to exhibit complete suppression 
of the electric dipole excitation in a narrow wavelength 
range when illuminated from the smaller disc side 11 . 
However, in a nanomaterial, there will be two counter- 
propagating waves and the electric dipole moment cannot 
be suppressed for both of them simultaneously 14 . The 
nanodimer geometry is described by the radii R\ = 15 nm 
and R2 = 20 nm of the discs and dimensions h = s = 
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10 nm defined in Fig. [5^. A cubic lattice with period 
A = 50 nm is now composed of the nanodimers such that 
the smaller discs are on the left-hand side. As a conse- 
quence of the asymmetry of the unit cell, we calculate 
the single layer response to illumination from both sides 
in order to obtain the reflection coefficients pL and pR. 

Using Eqs. (|5j) and (|9j) we calculate the transmission 
and reflection spectra for a nanodimer slab illuminated 
from the two sides and compare them with the numeri- 
cal results. Figure [5]o shows that while the theory very 
accurately resolves all spectral features, there is a slight 
deviation of the analytically obtained values from the ex- 
act numerical values for the wavelengths around 600 nm. 
This deviation obviously originates from the evanescent- 
wave coupling between the adjacent layers. However, 
considering that the gap size between the discs in the 
adjacent layers is only 20 nm, the agreement is still re- 
markably good. We obtained a similar good agreement 
between our theory and the numerical calculations also 
for angles of incidence of 0, 30 and 60 degrees as well as 
for the TE-polarization. The reflection spectra of 1000 
layers of nanodimers are depicted in Fig. [5^. When illu- 
minated from the small disc side, the reflection coefficient 
significantly decreases at around the electric dipole sup- 
pression wavelength of 618 nm. 

The accuracy of the presented theory depends on the 
extent of the evanescent waves in the direction of light 
propagation. Qualitatively, for the theory to be exact, 
the evanescent waves associated with the cut-off diffrac- 
tion orders must have a decay length that is much shorter 
than the spacing between the particles in two adjacent 
layers. For a two-dimensional square array of period 
A x = A yi the longitudinal wave vector of the first such 
order k z \ = [k 2 — (2tt/A x — k x ) 2 } 1 ^ 2 is imaginary, with 
k and k x being the magnitudes of the total and trans- 
verse wave numbers of the incident light. The next order 
would be k z2 = [k 2 - (2tt/A x + /^) 2 ] 1/2 - We define the 
decay length S of the evanescent field to be the distance 
for which the field amplitude has decayed by a factor 
of exp(— 1). The gap d between the particles in the z 
direction must then be much larger than 
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Considering the nanoshells with A^ = 130 nm, Ao = 
500 nm and k x = k/y/2, we obtain a decay length 
5 ~ 34 nm that is smaller than the 60 nm gap be- 
tween the adjacent shells. For the nanodimers, with 
A x = 50 nm, Ao = 500 nm and k x = k/y/2, the de- 
cay length is 5 ~ 9 nm. This value is smaller than 
the gap d = 20 nm between the nanodimers, which sup- 
ports the success of our analytical calculations. We note 
that for a nanomaterial to be homogenizable, it is re- 
quired that k « 2tt/A x . Equation ([TO]) then yields an 
evanescent wave decay length of A x /(2tt) and the crite- 
rion for neglecting the evanescent- wave coupling becomes 
d » A x / (2tt). As a practical criterion for when our the- 
ory can be applied to a homogenizable nanomaterial, we 
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FIG. 5. (a) Geometry of the disc nanodimer [Ri — 15 nm, 
R2 — 20 nm, h = s = 10 nm, A = 50 nm]. (b) is as in Fig. [3Jd. 
The green (red) color corresponds to the illumination from 
the side of the smaller (larger) discs, (c) Reflection spectra 
for 1000 layers as obtained from Eqs. ([5]) and ©. 



require that d > A x /2. 

In order to verify the above predictions on the influ- 
ence of the evanescent waves, we numerically calculate 
the transmission through the nanodimer slab, while vary- 
ing the transverse and longitudinal periods separately. In 
Fig. [6^i the transmission coefficient is plotted for an in- 
creasing longitudinal period A z . The transverse period 
is fixed to A^ = 100 nm in order to have the evanescent- 
wave coupling significant enough to cause a discrepancy 
at small A z between the theory and the numerical results. 
We notice that when d exceeds A x /2 [A z exceeds 80 nm], 
this discrepancy disappears. The transmission coefficient 
for an increasing transverse period A^ = A y is shown in 
Fig. [6}d. In this case the longitudinal period is fixed to 
A z = 70 nm, such that the evanescent-wave coupling is 
negligible at A x = 50 nm. As the transverse period is 
increased, the discrepancy between the theory and the 
numerical results starts to appear due to an increase in 
the evanescent wave decay length. We notice however 
that as A x is increased, the array also gets sparse, which 
reduces the influence of the nanoscatterers on the prop- 
agating wave. This effect counterbalances the growing 
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FIG. 6. Transmission coefficient t of a five- layer thick nan- 
odimer slab for (a) increasing period A z with fixed A x = A y = 
100 nm and (b) increasing A x — A y with fixed A z — 70 nm. 
A TM-polarized illumination with Ao = 600 nm and — 45° 
(from smaller disc side) is considered. The numerically calcu- 
lated results (solid lines) are shown along with the analytical 
results (stars) obtained from Eqs. ((5} and J9}. 



decay length of the evanescent waves such that the dis- 
crepancy between the theory and the numerical results 
remains small. 

In summary, we have introduced a simple analyti- 
cal theory for the description of light interaction with 
optical nanomaterials. Recognizing the subwavelength 
size of the nanomaterial's unit cells, we found that the 
evanescent- wave coupling between adjacent monolayers 



of the nanomaterial does not significantly influence the 
light propagation in the material. As a consequence, 
the nanomaterial can be accurately described in terms 
of the plane-wave transmission characteristics obtained 
for a single isolated monolayer. The presented numerical 
calculations demonstrate the wide applicability of this 
remarkably simple analytical model. 

In contrast to existing theoretical approaches, our one 
also correctly describes three-dimensional arrays of bifa- 
cial nanoscatterers, which is of practical importance for 
a large variety of nanomaterials, such as metamaterials 
with asymmetric unit cells. For both symmetric and bi- 
facial nanomaterial slabs, our method enables rapid ex- 
traction of the transmission and reflection coefficients, 
which can for example be used to calculate the effective 
wave parameters, such as e and /i, via the known re- 
trieval procedures*^. Furthermore, propagation of opti- 
cal beams through such nanomaterials can be described 
by using the angular spectrum representation with our 
model applied to each plane- wave component. We be- 
lieve that the presented theory has the necessary simplic- 
ity and accuracy to accelerate the development of optical 
nanomaterials tailored for real applications. 
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